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Abstract 

Using the recent experimental data oi B ^ D^*\Tr,p), B D^*^Di*\ B — > 
J/'ijjK^*^ and various model calculations on form factors, we re-analyze the 
effective coefficients ai and a2 and their ratio. QCD and electroweak penguin 
corrections to oi from B D^*^D^s'^ and 02 from B J /ipK^*^ are estimated. 
In addition to the model-dependent determination, the effective coefficient ai 
is also extracted in a model-independent way as the decay modes B D^*^h 
are related by factorization to the measured semileptonic distribution of i? — > 
D^*^£i> at = m\. Moreover, this enables us to extract model-independent 
heavy-to-heavy form factors, for example, FfP^(m^) = 0.66 it 0.06 it 0.05 and 
Aq^* {m^) = 0.56 lb 0.03 it 0.04. The determination of the magnitude of 02 
from B J/tpK^*'^ depends on the form factors Ff^, Af2* and V^^* at 
= in^j/^- By requiring that 02 be process insensitive ( i.e., the value of 
02 extracted from J/tpK and J/ipK* states should be similar), as implied by 
the factorization hypothesis, we find that B K^*^ form factors are severely 
constrained; they respect the relation Ff^ (m'^j^^^) ps 1.9Af ^* (rrij^^). Form 
factors A2^ and V^^ at q"^ = i^^'j/^ inferred from the measurements of the 
longitudinal polarization fraction and the P-wave component in 5 — > J /ipK* 
are obtained. A stringent upper limit on a2 is derived from the current bound 
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on B ^ D^TT and it is sensitive to final-state interactions. 



Typeset using REVI^ 



I. INTRODUCTION 



Nonleptonic two-body decays of B and D mesons have been conventionally studied in 
the generalized factorization approach in which the decay amplitudes are approximated by 
the factorized hadronic matrix elements multiplied by some universal, process-independent 
effective coefficients af^. Based on the generalized factorization assumption, one can catalog 
the decay processes into three classes. For class-I decays, the decay amplitudes, dominated 
by the color-allowed external H^-emission, are proportional to ai''^(Oi)fact where Oi is a 
charged current-charged current 4-quark operator. For class-II decays, the decay amplitudes, 
governed by the color-suppressed internal VF-emission, are described by a|''^(02)fact with O2 
being a neutral current-neutral current 4-quark operator. The decay amplitudes of the class- 
Ill decays involve a linear combination of a^*^(Oi)fact and a2^(02)fact- If factorization works, 
the effective coefficients a^^ in nonleptonic B 01 D decays should be channel by channel 
independent. Since the factorized hadronic matrix elements (Oi)fact are renormalization 
scheme and scale independent, so are af^ . 

What is the relation between the effective coefficients and the Wilson coefficients in 
the effective Hamiltonian approach ? Under the naive factorization hypothesis, one has 

1 1 
=ci (//) + — C2(/i), a2(//) = C2(//) + — ci(//), (1.1) 

for decay amplitudes induced by current-current operators Oi,2(At), where ci,2(a*) are the 
corresponding Wilson coefficients. However, this naive factorization approach encounters 
two principal difficulties: (i) the above coefficients are scale dependent, and (ii) it fails to 
describe the color-suppressed class-II decay modes. For example, the predicted decay rate 
of — s> by naive factorization is too small compared to experiment. Two different 

approaches have been advocated in the past for solving the aforementioned scale problem 
associated with the naive factorization approximation. In the first approach, one incorporates 
nonfactorizable effects into the effective coefficients [1-3]: 

af = ci(/x) + C2(/x) + xM) , «f = C2(/x) + ci(/x) + X2(/^)) , (1.2) 

where nonfactorizable terms are characterized by the parameters Xi- Considering the decay 
B^ D+TT" as an example, xi is given by 
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where 

(SD,7r)_ 



{D+\{ch),_,\B^){ii-\{du),_,\Q) 

{D+\{cb),_jB'){n-\{du),_,\0) ' 
are nonfactorizable terms originated from color signlet-singlet and octet-octet currents, re- 
spectively, {qiq2)v-A = ^i7/.(l-75)g2, and {qi\''q2)v-A = Qi^''lf^i'^-l5)q2- The fx dependence 
of the Wilson coefficients is assumed to be exactly compensated by that of Xiit^) [4]- That is, 
the correct n dependence of the matrix elements is restored by second approach, 

it is postulated that the hadronic matrix element {0{^)) is related to the tree- level one via 
the relation {0{^)) — g{ii){0)tree and that g{ii) is independent of the external hadron states. 
Explicitly, 

C(/X)(0(/X)) = C{pi)g{pi){0)tree = {0)tree. (1-5) 

Since the tree-level matrix element (O)tree is renormalization scheme and scale independent, 
so are the effective Wilson coefficients cf^ and the effective parameters af^ expressed by [5,6] 

af = cf + cf (i- + xi) , af = cf + cf (i- + X2) ■ (1-6) 

Although naive factorization does not work in general, we still have a new factorization 
scheme in which the decay amplitude is expressed in terms of factorized hadronic matrix 
elements multiplied by the universal effective parameters af^2 provided that xi,2 are univer- 
sal (i.e. process independent) in charm or bottom decays. Contrary to the naive one, the 
improved factorization scheme does incorporate nonfactorizable effects in a process indepen- 
dent form. For example, Xi = X2 = — | in the large- A^^^ approximation of factorization. 
Theoretically, it is clear from Eqs. (1.3) and (1.4) that a priori the nonfactorized terms Xi 
are not necessarily channel independent. In fact, phenomenological analyses of two-body 
decay data of D and B mesons indicate that while the generalized factorization hypothesis 
in general works reasonably well, the effective parameters af^2 do show some variation from 
channel to channel, especially for the weak decays of charmed mesons [1,7]. However, in the 



4 



energetic two-body B decays, Xi ^-re expected to be process insensitive as supported by data 
[4]. 

The purpose of the present paper is to provide an updated analysis of the effective 
coefficients af^ and af^ from various Cabibbo-allowed two-body decays of B mesons: B — > 
D^*)D^*\ D^*\7r,p), J/tpK'^*\ It is known that the parameter \af^\ can be extracted from 
L){*)+(7r-,p-) and B^ D(*W(*\ \af\ from B J/^jK^*), B^ L>W07rO(p°), 
and a^/a'f from B~ — > D'^*\'K,p). However, the determination of a^^ and is subject 
to many uncertainties: decay constants, form factors and their dependence, and the 
quark-mixing matrix element Vch- It is thus desirable to have an objective estimation of 
a^2- ^ model-independent extraction of ai is possible because the decay modes B D^*^h 
can be related by factorization to the measured semileptonic decays B — > D^*Hp. As a 
consequence, the ratio of nonleptonic to differential semileptonic decay rates measured at 

— ml is independent of above-mentioned uncertainties. The determination of |a|^| from 
B J/tpK^*^ is sensitive to the form factors Ff^, Af^' and V^^' at = mj/^. In 
order to accommodate the observed production ratio R = V{B J/ipK*)/T{B — > J/ijjK) 
by generalized factorization, af^ should be process insensitive; that is, af^ extracted from 
J/ipK and J/ipK* final states should be very similar. This puts a severe constraint on the 
form-factor models and only a few models can satisfactorily explain the production ratio R. 

The rest of this paper is organized as follows. In Sec. II, we introduce the basic formula 
and the classification of the relevant decay modes which have been measured experimentally. 
Sec. Ill briefly describes various form-factor models. The results and discussions for the 
effective parameters a'^ and are presented in Sees. IV and V, respectively. Finally, the 
conclusion is given in Sec. VI. 



II. THE BASIC FRAMEWORK 



Since, as we shall see below, the decays B — > D^*^D^*\ J/'^K^*'> receive penguin contri- 
butions, the relevant AS = 1 effective Hamiltonian for our purposes has the form 



Kff - ^|K,v;;[ci(/.)or(/.) + c2(/.)or(/.) 

10 ^ 
i=3 J 



ci(/.)oS-)(/.)+c2(/.)or(/^)j 

(2.1) 



where 



0^^ = (c6),_,(g«),_,, 

03(5) = (#)v_aE(^V)v_a(v+a)> 



(s6)^_^(cc)^_^, 



04(6) = {(lah0)y_j, J2Wa) 



v-a{v+a}i 



(2.2) 



08(10) = 7:{<iabp)y_^Y.^i'Wo 



a)v+A{v-A)i 



with O3-O6 being the QCD penguin operators and O7-O10 the electroweak penguin opera- 
tors. 

To evaluate the decay amphtudes for the processes B — > D^*^D^*\ D'^*)+'°(7r^, p~), 
J/ipK^*\ we first apply Eq. (1.5) to the effective Hamiltonian (2.1) so that the factorization 
approximation can be applied to the tree-level hadronic matrix elements. We also introduce 
the shorthand notation X^^^^'^^^ to denote the factorized matrix element with the F2 meson 
being factored out [6], for instance, 



(2.3) 



The results are: 



Class I: 7r-(p-) 

The decay amplitudes are given by 



^(s',^DW+7r-(p-)) 



Gp 
V2 



j^{B°DM+,n-{p-)) ^ ^^^{B°,DM+n-{p-)) 



(2.4) 



where x^^'-^^*''^^ is the factorized H^-exchange contribution. 



Class I: B- D(*)Od(*)" and ^ D^*^- 
The decay amplitudes are given by 



A{B ^DD,) = ^|v;6K: «i - VtbVZ 



(24 + Ctio 



2(a6 og) 



m 



{mb-mc}{mc + ms 



(2.5) 



where use of V^tT^* = —VcbV*g has been made and 



ai[B ^ D Ds) = Oi 1 H h 2 



(2.6) 



Likewise, 



(2.7) 



Note that the decay B' D^D^/^' also receives a contribution from the W- 
annihilation diagram, which is quark-mixing-angle doubly suppressed, however. 

Class II: 5° ^ 7rO(pO) 

The factorized decay amplitudes are given by 



A{B: i5WV(p°)) = -^V,,V:, a. 



■^(B°,rO(pO),D(*)0) ^ ^(B°,D(*)07rO(pO)) 



(2.8) 



where X*^^ ,o(*'°7r°(p°)) jg ^j^g factorized VT-exchange contribution. 



• Class II: B+ J/i/j K^*^+ and S° ^ J/V' K^*^^ 
The decay amplitudes are given by 



v2 



where 



02(5 ^ J/^/T^*)) = 02 



cts + 05 07 ag 



02 



(2.9) 



(2.10) 



Class III: B- t^~{p~) 
The decay amplitudes are given by 



A(S-^i^W° 7r-(p-)) = ^V;,V:, 



ai 



^(B-DW0,7r-(p-)) ^^^^(B-7r-(p-),D(*)0)l _ ^^2.11) 



Class III: B- D^^^'^R- 

The factorized decay amplitudes are given by 



ai 



(2.12) 
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Under the naive factorization approximation, a2i — c|f + -^c^^i and 021-1 = c^-i + 

{i — 1, • • • , 5). Since nonfactorizable effects can be absorbed into the parameters af^, this 

amounts to replacing Nc in Oj by {Nf^)i [6] with 

^^+Xi- (2.13) 



c 



{Nf)i K 
Exphcitly, 

„eff _ eff , ^ ^eff eff _ eff , ^ ^eff 14^ 

(For simphcity, we have already dropped the superscript "efF' of in Eqs. (2.4) to (2.12) 
and henceforth.) 

Although the purpose of the present paper is to treat the effective coefficients Oi and 02 
as free parameters to be extracted from experiment, it is clear from Eqs. (2.6) and (2.10) that 
the determination of ai,2 from B — > £)(*)£)(*) and J/t/jK^*^ is contaminated by the penguin 
effects. Therefore, it is necessary to make a theoretical estimate on the penguin contribution. 
To do this, we employ the effective renormalization-scheme and -scale independent Wilson 
coefficients c^^ obtained at k'^ = ml/2 {k being the gluon's virtual momentum) [6]: 

cf = 1.149, cf = -0.325, 

cf = 0.0211 + ^0.0045, cf = -0.0450 - i0.0136, 

cf = 0.0134 + i0.0045, cf = -0.0560 - i0.0136, 

cf = -(0.0276 + i0.0369)Q;, cf = 0.054 a, 

cf = -(1.318 + i0.0369)a, cf^ = 0.263 a. (2.15) 

For nonfactorizable effects, we choose N'f^LL) ^ 2 (see Sec. V.E) for {V — A){y — A) 
interactions (i.e. operators Oi, 2.3,4,9,10) and N^^{LR) ^ 5 for {V — A){y + A) interactions 
(i.e. operators 05^6,7,8)- Our choice for Nf^{LR) is motivated by the penguin-dominated 
charmless hadronic B decays (for details, see [6,8]). Hence, the theoretical values of the 
effective coefficients are given by 

ai = 0.986, a2 = 0.25, 

03 = -(0.00139 + 0.00226i), 04 = -(0.0344 + 0.0113i), 

as = 0.0022 + 0.00181i, Og = -(0.0533 + 0.0127i), 



a7 = -(1.24 + 2.73i) x lO'^ 
ag = -(87.9 + 2.73i) x 10'^, 



as = (3.59 - 0.55i) x 10-^ 
= -(29.3 + 1.37i) x 10"^ 



(2.16) 



Prom Eqs. (2.6), (2.7), (2.9) and (2.16), penguin corrections to the tree amplitudes are found 
to be ^ 



\Ap/At\{B^DD,)^ 

\Ap/At\{B^D*D,)^ 
\Ap/At\{B dWd;) = 
\Ap/AT\{B^J/tl;K^*^)^ 



m 



Ds 



ai 

(24 + aio 



ai {rrib - mc)ims + mc) 
ae + as 



Oi Oi (rrib + mc)(ms + rric) 

04 + Ctio 



= 0.159, 
= 0.037, 



0.040, 
03 + 05 + 07 + Og 



02 



0.033 , 



(2.17) 



where we have used the current quark masses at the scale /i — 0(m6): mslrrih) — 
105 MeV, mc{mb) — 0.95 GeV, mb{mb) — 4.34 GeV. Therefore, the penguin contribution 
to B ^ D*Ds, D^*W* and J/i/jK^*^ is small, but its effect on B ^ D Ds is significant. 
Numerically, the effective di defined in Eqs. (2.6) and (2.10) are related to by 

di{B^DDs) = 0.847 ai 
61 (S ^ D*Ds) = 1.037 Oi 
di{B D^*^D;) = 0.962 oi 
02(5 ^ J/jjjK^*^) = 0.968 02. (2.18) 

To evaluate the hadronic matrix elements, we apply the following parametrization for 
decay constants and form factors [10] 



{0\A^\P{q)) = tfpq^,, {0\V,\V{p,e)) = fymvs^, 



{P'{p')\V,\P{p))^i^p,+p'^ 
{V{p',e)\V,\P{p)) ^ 



m p — ml 



2>. TTlp Tflpt 



qn, 



nip + iriv 



e,.ape*W^V{q'), 



{V{p',e)\A^\P{p))^i\{mp + mv)e^Ai{q'') - 



e ■ p 
mp + mv 



{p + p'),A,{q') 



^Our numerical estimate for the penguin effects m. B ^ D Dg differs from [9] due to different 
choices of N^^{LL), N^^{LR) and running quark masses. 
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£ ■ p 



-2mv qM^ir) - Mr)] 



where q^p-p', Fi(0) = Fo(0), ^3(0) = A(0), 

Asiq') 



mp + my 2n mp-mv . . 2-. 

A(g ) ^ A2{q ) , 



(2.19) 



(2.20) 



2mv 2mv 
and P, V denote the pseudoscalar and vector mesons, respectively. The factorized terms in 
(2.4)- (2. 12) then have the expressions: 

Xm,P.) ^ (P2|(g2g3)^_jo)(Pi|(gi6),_jE) = tfp.iml - mljFo^^HO, 
X(^^y^ = {V\iq2q3)y_jO){Pm),_jB) = 2fymvFr{m'y){e*-p,), 
Xi^y^n ^ {P\(q,q,)^_JO){V\{q,b),_jB) = 2fpmyA^''iml){e*.p,), 



XiBV^V.) ^ {V2\{q2qs)y_JO)mq,b),_jB) = -tfy. 



2A^^Hm: 



[mB + mij 



(2.21) 



{niB + mi) 

where e* is the polarization vector of the vector meson V. 

With the factorized decay amplitudes given in Eqs. (2. 4) -(2. 12), the decay rates for 
B — > PP, VP are given by 

Pc 



r{B ^ P1P2) 



pI 



\A{B^PiP2)\\ 



T{B VP) = ^\A{B VP)/{e.p,)\\ 



(2.22) 



where 



Pc 



\J[rrrB - [nii + /;;2)-][/"i' - ("'1 - "''2)' 
2mB 



(2.23) 



is the cm. momentum of the decay particles. For simplicity, we consider a single factorizable 
amplitude for B VV: A{B V1V2) = aX^^^i-^^) xhen 



r{B ^ V1V2) 



Pc 



87rm| 



\a 



'mB + mi)m2fv,Af'''{ml)\''H, 



with 



H = (a-bx)' + 2{l + cY), 



(2.24) 



(2.25) 



and 
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2m%pl 



2mBPc 



2mim2 



X 



mim2{mB + mi)^ ' 



(niB + 



[ml) 



(2.26) 



where mi (777.2) is the mass of the vector meson Vi (V2). 

Prom Eqs. (2.4-2.11) we see that |ai| can be determined from ^ D^*'>~^{n~, p~), B 
D^*^D^;\ lasi from B J /ijK^*\ 5° L)(*)°(7r°, p°), provided that the VT-exchange 
contribution is ncghgible in B ^ D'^*\7r,p) decays and that penguin corrections are taken 
into account. It is also clear that the ratio 02/01 can be determined from the ratios of charged 
to neutral branching fractions: 



i?i = 



B{B- 



71 



R9. 



B{B" 
B{B- 



D+TT- 



Ra = 



~ i3(5° 
B{B~ - 



► D+p- 



i3(B° 
B{B- 



D*+n- 
> D*^p- 



B{W D*+p- 



t(B- 



r{B-) 

r{B-) 
r(B') 

rjB-) 



1 + 



1 + 



1 + 



m 



B 



mlfp F^^-jml) a2 

f,Ff^{ml)a^) ' 
U Ar\ml)aJ ' 



Hi 



H 



H 



with 



Hi 
H2 



m 



D* 



a2 



m 



,{mB + mD*) fp Af^*(m2) oi' 
(a - bx){a - h'x') + 2(1 + cc'yy'), 

(a-6V)2 + 2(l + cV), 



(2.27) 



(2.28) 



where H,a,b,c,x,y are those defined in Eqs. (2.25) and (2.26) with Vi = D* and V2 = p, 
and b\c\x\y' are obtained from 6, c, a;, y respectively with the replacement D* p. 



III. MODEL CALCULATIONS OF FORM FACTORS 



The analyses of 02, oi, and 02/01 depend strongly on the form factors chosen for cal- 
culations. In the following study, we will consider six distinct form-factor models: the 
Bauer-Stech-Wirbel (BSW) model [10,11], the modified BSW model (referred to as the 
NRSX model) [12], the relativistic light-front (LF) quark model [13], the Neubert-Stcch 
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(NS) model [4], the QCD sum rule calculation by Yang [14], and the hght-cone sum rule 
(LCSR) analysis [15]. 

Form factors in the BSW model are calculated at zero momentum transfer in terms 
of relativistic bound-state wave functions obtained in the rclativistic harmonic oscillator 
potential model [10]. The form factors at other values of are obtained from that at = 
via the pole dominance ansatz 

where rUpoie is the appropriate pole mass. The BSW model assumes a monopole behavior 
(i.e. n = 1) for all the form factors. However, this is not consistent with heavy quark 
symmetry for heavy-to- heavy transition. In the heavy quark limit, the B ^ D and B D* 
form factors are all related to a single Isgur-Wise function through the relations 



I^Jm^nvE, 1 - q^/{mB + 1710^ 



— . (3.2) 

1 - q"^ / {tub + mD*y 

Therefore, the form factors Fi, V, ^0,^2 in the infinite quark mass hmit have the same 
dependence and they differ from Fq and Ai by an additional pole factor. In general, the 
heavy-to-heavy form factors can be parametrized as 

ABD*j2-._ ( + rnp* .BD'j(X. ^B - ^D* aBD' fr.^\ ^' il ) 

aBD*(^2^ ( ^3 + mp* \ ujp* jq^) + 1 ^^^2^ 

where 
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2mBmj)(,) 



G{Q')-G{Qla.)[^-pU^D{q')-l)], 
Hq')^ ^(?LJ[1 - Pd*Kw(?') - 1)] , 



r{q^)-- 
ri{q^)-- 
r2{q^)-- 



^'{<iL.)[^-p'\^\<i')-^)] 



(niB + moY 



{niB + mD*Y 
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Fo^^iq') ' 



(3.4) 



{niB + ruD'Y _ 

In the heavy quark hmit oo, the two form factors J^{q^) and G{q^)i whose slopes are 

coincide with the Isgur-Wise function ^{q^)^ and r(g^), ri(g^) as well as r2{q^) are equal 
to unity. The q^ dependence of S ^ D^*-* form factors in the NRSX and NS models is more 
complicated because perturbative hard gluon and nonperturbative l/mg corrections to each 
form factor are taken into consideration and moreover these corrections by themselves are 
also q^ dependent (see [12] for more details). 

Form factors for heavy-to-heavy and heavy-to-light transitions at time-like momentum 
transfer are explicitly calculated in the LF model. It is found in [13] that the form factors 
Fi,V,Aq,A2 all exhibit a dipole behavior, while Fq and Ai show a monopole dependence 
in the close vicinity of maximum recoil (i.e. (f' — 0) for heavy-to-light transitions and in 
a broader kinematic region for heavy-to-heavy decays. Therefore, the dependence of 
B — > form factors in the heavy quark limit is consistent with the requirement of heavy 
quark symmetry. Note that the pole mass in this model obtained by fitting the calculated 
form factors to Eq. (3.1) is slightly different from that used in the BSW model (see Table I). 

Due to the lack of analogous heavy quark symmetry, the calculation of heavy-to-light 
transitions is rather model dependent. In addition to the above-mentioned BSW and LF 
models, form factors for the B meson to a light meson are also considered in many other 
models. The NRSX model takes the BSW model results for the form factors at zero momen- 
tum transfer but makes a different ansatz for their q^ dependence, namely a dipole behavior 
(i.e. n = 2) is assumed for the form factors Fi, Aq, A2, V ^ motivated by the heavy-quark- 
symmetry relations (3.2), and a monopole dependence for Fq^Ai. The heavy-to-light form 
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factors in the NS model have the expressions [4] : 



fBP(2^ ( mB + mp \ \/ ^Bp(g^)+T 



Arid') 



<^sp(0)-l 

itlb + iTLp I 2 1 u;bp(0) — C(;bp(^^- 



2^mBmp\l ujBpiq^) + 1 1 + upp^q^) - ujBp{jn\^) 
mp + niv I 2 1 cc;bv(0) - cjBV^("^io-) 



1^2^/^7I^yl V 2 I + rv '^Bvi^)-! ' 



mB + niv j 2 1 cjBy(O) - cjbv("^i+) 



where 



^y- — nr-z — 1 1 + ^2 _^^„_^.^2 • (3-6) 



2msmBP(y) 
(me — my)^ / ^ AniBmv 
AmBmv \ m\+ — {niB — mvY 

Here mo-, mi-, and mi+ are the lowest resonance states with the quantum numbers 0~, 1~, 

and 1+, respectively.^ 

We consider two QCD sum rule calculations for S-to-hght transitions. The form factors 
Fi and Ai in the Yang's sum rule have a monopole behavior, while A2 and V show a dipole 
dependence. The momentum dependence of the form factors p^'^^^^ and Aq"^^^ ^ is slightly 
complicated and is given by [14] 

F„'"<(g^)= -0.32 ( -r^] + F^if) , 



^5.8^ — q^ ) rriB — ttik 



^ 1 5.982 -g2/ I 2m. ^ ^ 2mp ^ 



(,^)= 0.02,^ (^ll^)' + (^^^l^-^rto^) - ^^^l^-^r-tf ) ) . (3.7) 



^Following [4], we will simply add 400 MeV to mi- to obtain the masses of 1+ resonances. 
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The behavior of S-to-hght form factors in the LCSR analysis of [15] are parametrized as 

^ _ , (3.8) 

where the relevant fitted parameters ap and hp can be found in [15]. 

Since only the form factors for S-to-light transition are evaluated in the Yang's sum rule 
analysis and the LCSR, we shall adopt the parametrization (3.3) for the B — > D^*'> form 
factors, in which the relevant parameters are chosen in such a way that B — > D^*^ transitions 
in the NS model are reproduced: 

^(gLj = 0.88, ^(gLj = 1.00, 

p^ = 0.62, p'2 = 0.62, p^. =0.91, 

r{q^)^l, ri(?^) ri = 1.3±0.1, r2(g^) r2 = 0.8 ± 0.2 , (3.9) 

as a supplement to the Yang's [14] and LCSR [15] calculations. The theoretical prediction 
for ri and r2 [16] is in good agreement with the CLEO measurement [17]: ri = 1.18 ± 0.32 
and r2 = 0.71 ± 0.23 obtained at zero recoil. Note that the predictions of S — > D^*'' form 
factors are slightly different in the NRSX and NS models (see Table I) presumably due to 
the use of different Isgur-Wise functions. 

To close this section, all the form factors relevant to the present paper at zero momentum 
transfer in various models and the pole masses available in the BSW and LF models and in 
the Yang's sum rules are summarized in Table I. 



IV. DETERMINATION OF ai 

In order to extract the effective coefficient ai from S° (tt", p~) and B D^*^D^*^ 

decays, it is necessary to make several assumptions: (i) the W^-exchange contribution in 
is negligible, (ii) penguin corrections can be reliably estimated, and (iii) 
final-state interactions can be neglected. It is known that H^-exchange is subject to helicity 
and color suppression, and the helicity mismatch is expected to be more effective in B 
decays because of the large mass of the B meson. Final-state interactions for B — > D^*\7r, p) 
decays are customarily parametrized in terms of isospin phase shifts for isospin amplitudes. 
Intuitively, the phase shift difference 5i/2 — ^3/2, which is of order 90° for D —>■ Kn modes, 
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is expected to play a much minor role in the energetic B — > Dn decay, the counterpart of 
D — > Ktt in the B system, as the decay particles are moving fast, not allowing adequate 
time for final-state interactions. Prom the current CLEO limit on B^ — > D°7r° [18], we find 



sin-(^v.-M<^ ;(^;B(g-^g„,-) =0.109. (4.1) 

and hence 

\Sl/2-Ss/2\B^Dn<l9°. (4.2) 

We shall see in Sec. V.III and in Fig. 1 that the effect of final-state interactions ^ subject to 
the above phase-shift constraint is negligible on r(B^ — > D+7r~) and hence it is justified to 
neglect final-state interactions for determining ai. The extraction of 02 from B J/ipK^*^ 
does not suffer from the above ambiguities (i) and (iii). First, H^-exchange does not con- 
tribute to this decay mode. Second, the J/ipK^*^ channel is a single isospin state. 



A. Model-dependent extraction 

We will first extract ai from the data in a model-dependent manner and then come to 
an essentially model-independent method for determining the same parameter. 

Armed with the form factors evaluated in various models for B D and B — >■ D* 
transitions, wc arc ready to determine the effective coefficient ai from the data of 5^ — > 
£)(*)+ (7r~,p") and B £)(*)!:)(*) decays [20]. The results are shown in Tables II and III in 
which we have taken into account penguin corrections to oi [see Eq. (2.18)]. We will choose 
the sign convention in such a way that Oi is positive; theoretically, it is expected that the sign 



^Final-state interactions usually vary from channel to channel. For example, \Si/2 — ^3/2! is of 
order 90° for D Ktt, if*7r, but it is consistent with zero isospin phase shift for D Kp. The 
preliminary CLEO studies of the heUcity amplitudes for the decays D*'^p~ and B~ D*^ p" 

indicate some non-trivial phases which could be due to FSI [19]. At any rate, FSl are expected to 
be important for the determination of the effective coefficient 02 (see Sec. V.III), but not for ai. 
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of oi is the same as ci. In the numerical analysis, we adopt the following parameters, quark- 
mixing matrix elements: \Vcb\ — 0.039±0.003, \Vud\ — \Vcs\ — 0.975±0.001; decay constants: 
/^=132 MeV, fK - 160 MeV, /p=216 MeV, /i?=200 MeV, /d*=230 MeV, /d.=240 MeV, 
/dj=275 McV, /j/^=394 MeV, and lifetimes: r(fi") = (1.57 ±0.03) ps, t{B-) = (1.67 ± 
0.03) ps [21]. Because of the uncertainties associated with the decay constants f^^ and fn*, 
the value of ai obtained from B —>■ D^*^D^*^ decays in Table III is normalized at fo^ — 240 
MeV and fo* — 275 MeV. For example, oi determined in the NRSX model reads 

ai D^*'>+ {n- , p-)) = 1.04 ±0.03 ±0.08, 

240 MeV\ 

ai (b ^ L»Wl>:) = (1.12 ± 0.12 ± 0.08) x j ^ (4.3) 

where the first error comes from the experimental branching ratios and the second one from 
the B meson lifetimes and quark-mixing matrix elements. Evidently, ai lies in the vicinity 
of unity. 

Several remarks arc in order, (i) From Tables II and III wc sec that ai extracted from 
B £)(*)£)(*) is consistent with that determined from B D^*^7r(p), though its central value 
is shghtly larger in the former, (ii) Theoretically, it is expected that r(B^ £)(*)+£)(*)-) — 
r{B- d(*)od(*)-) and hence B{B- d(*)od(*)-) ^ 1.07B(S° ^ dW+D(*)-). The 
errors of the present data are too large to test this prediction, (iii) The central value of Oi 
extracted from B^ D*^Dj and B" D*^Dj in the BSW model deviates substantially 
from unity. This can be understood as follows. The decay amplitude of the above two modes 
is governed by the form factor Aq^* {m'jjj. However, the dependence of ^o(?^) in this 
model is of the monopole form so that does not increase with fast enough compared 
to the other form-factor models. 

B. Model-independent or model-insensitive extraction 

As first pointed out by Bjorken [22], the decay rates of class-I modes can be related under 
the factorization hypothesis to the differential semileptonic decay widths at the appropriate 
q^. More precisely. 



ai 



;i.26± 0.11 ±0.09) 



X 
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Ct(*) _ 



(*) 



(4.4) 



where ai = ai in the absence of penguin corrections [the expressions of Oi are given in 
Eq. (2.6)], Vij = for /i = tt, p, Vij = for h = and [12] 



Yp 



ip — 
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J:^=o,±l\Hr*in^l)\ 
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Yv = Y*= 1 



(4.5) 



with the hehcity amphtudes i?o(?^) and H±{q^) given by 



Hl^\q') = {mB + mD*)Af^\q') T 



BD* , 



2mBPc 



1 



2mD*Vq^ 



[m. 



m 



D* 



TUB + mo* 

-q'){mB + mD^)Af''\q') 



(4.6) 



rriB + mo* 
where Pc is the cm. momentum. 

Since the ratio 5"^*^ is independent of Vcb and form factors, its experimental measurement 
can be utihzed to fix ai in a model-independent manner, provided that Yji ' is also indepen- 
dent of form-factor models. From Table IV we see that Y* and in particular Y.^ are essentially 
model independent. The BSW model has a larger value for Yd^ and a smaller value for Y^^ 
compared to the other models because all the form factors in the former are assumed to have 
the same monopole behavior, a hypothesis not in accordance with heavy quark symmetry. 
In the heavy quark limit, one has Yd^ 1.36 and Yg^ fa 0.37 [12]; the former is quite close 
to the model calculations (see Table IV). In short, Y^*^ ^ Y^} , Y„ are model independent, Y* 
is model insensitive, while Yd^ and Y^^ show a slight model dependence. 

In Table V the experimental data of dB{^ D^i^i?)/dq'^ (at q^ = and and 
dl3(^ D*+e-u)/dq^ are taken from [23] and [24], respectively. Note that the "data" 
of dB{B^ D'^l~V)/dq^ at small q^ are actually obtained by first performing a fit to the 
experimental differential (f distribution and then interpolating it to q^ — m^ and m^. For 
the data of dB(B^ — > D'^£~P)/dq^ at q^ — mj^^ and mjj* we shall use the CLEO data for 
dV/duj expressed in the form [25] 
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"^^^^Z^^^^ = ^ (^B + mnrmUu^' - 1)^/^1 K,^(a;)p, (4.7) 

where cu = vb ■ vd — {fn^B + "^l? ~ ?^)/(2mBmi3). A fit of J-{uj) parametrized in the hnear 
form 

r{u)^T{l)[l-p\uj-l)l (4.8) 

to the CLEO data yields [25] 

= 0.81 ± 0.14, \Vcb^{l)\ = (4.31 ± 0.42) x 10"^ (4.9) 

From (4.7)-(4.9) we obtain dB(K D+i-iy)/dq^ at = mjy and mj^* as shown in Table 
V. Note that we have applied the relation B{B~ D^*^^Di*'^-) ^ 1.07 S(5° D^*^+Di*^^) 
to get the average branching ratio for B £)(*)!)(*) and the ratios S^^] and S*^]. It is 
easy to check that the data, say dB/dq^ — (0.35 ± 0.06) x 10~^ GeV~^ at — ml, are well 
reproduced through this interpolation. 

The results of ai extracted in this model-independent or model-insensitive way are ex- 
hibited in Table V (for a recent similar work, see [26]), where we have chosen Yo^ = 1.36 and 
Y^^ = 0.40 as representative values. As before, the value of ai obtained from B [)(*)[)(*) 
decays is normalized at fo^ = 240 MeV and fo* = 275 MeV. In view of the present theoret- 
ical and experimental uncertainties with the decay constants fo^ and fo* and the relatively 
small errors with the data of Dn and D*7r final states, we beheve that the results (see Table 
V) 

ai(s" ^ D+TT") = 0.93 ±0.10, ai(# ^ D^+tt") = 1.09 ± 0.07 (4.10) 

are most reliable and trustworthy. Of course, if the factorization hypothesis is exact, Oi 
should be universal and process independent. However, we have to await more precise 
measurement of the differential distribution in order to improve the values of Oi and to have 
a stringent test on factorization. 

Once Oi is extracted from Sf^\ some of the B — > D^*^ form factors can be determined 
from the measured B — > D^*\'K,p) and £)(*)£)(*) rates in a model- independent way: 

F^^{ml) = 0.66 ± 0.06 ± 0.05 , ^o^^(^dJ = 0.78 ± 0.08 ± 0.06 , 

Ff^(ml) = 0.67 ±0.06 ±0.05, Ff^(ml.) = 0.89 ± 0.10 ± 0.07 , 

A^^*(m^) = 0.56 ±0.03 ±0.04, A^^* (m^J = 0.77 ± 0.03 ± 0.06 . (4.11) 
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It should be stressed that the above form-factor extraction is independent of the decay 
constants fo^ ^-nd fo* ■ It is interesting to see that Ff^ tends to increase with q'^ faster than 
F^^, in agreement with the heavy-quark-symmetry requirement (3.2). 

The decay constants fu^ and fu* can be extracted if ai determined from B D^*^D^/'> 
is assumed to be the same as that from D^*^7r(p) channels. For example, the assumption of 
ai(B^ D^D~) — ai{B^ D~^tt~) will lead to an essentially model-independent deteTmi- 
nation of fo^ ■ We see from Table V that 

(1.01 ±0.14) (240 MeV/ZijJ = 0.93±0.10, (4.12) 

and hence 

fo, = (261 ± 46) MeV. (4.13) 

Another equivalent way of fixing fo^ is to consider the ratio of hadronic decay rates [4] 

g(# ^ D^D;) ^ ( FrimlJ fnA ' 1.812 GeV 

B(B'^D+n-) [ F,^^{ml) fj 2.306GeV' ^^"'^^ 

where 1.812 GeV and 2.306 GeV are the cm. momenta of the decay particles Dg and tt, 
respectively, use of ai{B — > DDg) — ai{B Dtt) has been made and penguin corrections 
have been included. It is easy to check that the same value of is obtained when the 
model-independent form factors (4.11) are applied to (4.14). Likewise, 

/d* = (266 ± 62) MeV (4.15) 

is obtained by demanding ai(B^ D'^D*~) — ai(B^ — > D'^p"), for example. However, it is 
worth stressing again that the above extraction of /d^ and fn* suffers from the uncertainty 
of using the same values of Oi for different channels [26]. Since the energy released to the 
DDg state is smaller than that to the Dtt state, Oi may differ significantly in these two decay 
modes. 

V. DETERMINATION OF 02 AND aa/ai 

In principle, the magnitude of 02 can be extracted directly from the decays B — > J/t/jK'^*^ 
and B^ D(*V{p^) and indirectly from the data of 5" ^ D^*h{p) and B^ D^*h{p). 
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Unfortunately, the branching ratios of the (class-II) color-suppressed decay modes of the 
neutral B meson are not yet measured. Besides the form factors, the extraction of 02 from 
B D^*^7r{p) depends on the unknown decay constants and fn*- On the contrary, 
the decay constant /j/,^ is well determined and the quality of the data for B J/ipK^*'' 
is significantly improved over past years. Nevertheless, the relative sign of oi and 02 can 
be fixed by the measured ratios • • • , i?4 [cf. Eq. (2.27)] of charged to neutral branching 
fractions of S — > D^*^n{p), and an upper bound on \a2\ can be derived from the current limit 
on S° ^ L>V. 



A. Extraction of \a2\ from B — > J/ijjK^*^ 

From Eqs. (2.9) and (2.21), it is clear that 02 derived from B J/ipK and B J/ipK* 
depends on the form factors 

pBK 

{rn?jl^) and A^^* (mjy^), V^^* (w?ji^). These form factors 
evaluated in various models are collected in Table VI. A fit of Eq. (2.9) to the data of 
B{B^ J/i/jK) (see Table VII) yields 

\a,\{B J/^K) = (0.26 ± 0.02) f ) • (5.1) 

Prom Table VII we also see that the extracted value of |a2|(-B — > J/t/jK*) in various models 
can be approximated by 

\a2\(B ^ J/V^X*) « (0.21 ± 0.02) ( ,J.f ^ J • (5-2) 

This implies that the quantity \/H defined in Eq. (2.25) is essentially model-independent, 
which can be checked explicitly. If the factorization approximation is good, the value of 
a2 obtained from J/ipK and J/ipK* states should be close to each other. This is justified 
because the energy release in S — > J/t/jK* is similar to that in S — > J/t/jK and hence the 
nonfactorizable effects in these two processes should be similar. However, we learn from 
Table VII that only the NRSX, LP models and the Yang's sum rule analysis meet this 
expectation. 

In order to have a process-insensitive a2, it follows from Eqs. (5.1) and (5.2) that the 
form factors Ff^ and Af^* must satisfy the relation 
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It is evident from Table VI that the ratio Ff^{m'^j^,^)/Ai^*{m'^j^^) is close to 1.9 in the 
aforementioned three models. This is also reflected in the production ratio 

Based on the factorization approach, the predictions of R in various form-factor models are 
shown in Table VIII. The BSW, NS and LCSR models in their present forms are ruled out 
since they predict a too large production ratio. To get a further insight, we consider a ratio 
defined by 

which measures the enhancement of Ff^/Af^* from — to finite q^. Z is close to unity 
in the BSW model and in Yang's sum rules (see Table VI) because Ff^ and Af^* there 
have the same monopole dependence, while in the other models 

pBK 

increases with q^ 

faster than Ai. For example, the dependence of F^^* in the LF model differs from that 
of Af^* by an additional pole factor. We see from Table VI that NS, LCSR and LF models 
all have similar q'^ behavior ^ for Z with Z(mj/^) ~ 0(1.35). In order to accommodate the 
data, we need Ff^{0)/Af^*{0) > L30. However, the values of Ff^(0) and >lf^*(0) are the 
same in both NS and LCSR models (see Table I) and this explains why they fail to explain 
the production ratio. By contrast, although Z 1 in the Yang's sum rules, its Ff^[0) is 
two times as large as Af^*(0) so that Ff^{m]/^)/Af^*{m'^j/^) ^ Ff^{0)/Af^*{0) ^ 2. 
We thus conclude that the data oi B —> J/ipK^*'> together with the factorization hypothesis 
imply some severe constraints on the B — > K^*^ transition: the form factor Ff^ must be 
larger than Af^* by at least 30% at g^ = and it must grow with q"^ faster than the latter 
so that Ff^{m'j/^)/Af^\m]/^) « L9. 

Since experimental studies on the the fraction of longitudinal polarization F^/F and the 
parity-odd P-wave component or transverse polarization |Pp measured in the transversity 



^Although Ff^ has the same dipole q'^ behavior in NRSX and LF models, its growth with q"^ in 
the former model is slightly faster than the latter because of the smaller pole mass. 
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basis in S — > J/ipK* decays are available, we have analyzed them in various models as shown 
in Table VIII. At this point, it is worth emphasizing that the generalized factorization hy- 
pothesis is a strong assumption for the B VV decay mode as its general decay amplitude 
consists of three independent Lorentz structures, corresponding to P and D-waves or 
the form factors Ai, V and A2. A priori, there is no reason to expect that nonfactorizable 
terms weight in the same way to S-, P- and D-waves. The generalized factorization as- 
sumption forces all the nonfactorizable terms to be the same and channel-independent [27]. 
Consequently, nonfactorizable effects in the hadronic matrix elements can be lumped into the 
effective coefficients Oj under the generalized factorization approximation. Since the decay 
B — > J/ipK^*^ is color suppressed and since |ci/c2| ^ 1, it is evident from Eq. (1.6) that 
even a small amount of nonfactorized term X2 will have a significant impact on its decay 
rate. However, it is easily seen that nonfactorizable effects are canceled out in the produc- 
tion ratio, the longitudinal polarization fraction and the P-wave component. Therefore, the 
predictions of these three quantities are the same in the generalized and naive factorization 
approaches. Explicitly [28,27], 



«=,08|. I^ = ^M. |Pr = Hf!. (5.6) 



where H,a,b,c,x,y are defined in Eqs. (2.25) and (2.26). Numerically, a = 3.165, b = 
1.308, c = 0.436. Form factors Af^* and V^^* at — m?ji^ can be inferred from the 
measurements of Fi/F and |Pp in B J/j/jK*. For illustration we take the central values 
of the CLEO data [29] (see also Table VIII): R = 1.45, F^/F = 0.52 and |Pp = 0.16. Since 
z 1.9, it follows from Eq. (5.6) that 



A^^'lm^u.) V 



IIL ■ 



x^ 7^^,; 1.19, \bk' r[ - 1-45 ■ (5.7) 

Prom Table VIII we see that all the model predictions for F^/F and |P|^ are in agreement 
with experiment ^ except that the longitudinal polarization fraction obtained in the NRSX 
model is slightly small. Indeed, among the six form-factor models under consideration, the 



^Historically, it has been shown [28] that the earlier data of R and Tl/V cannot be simultaneously 
accounted for by all commonly used models for form factors. In particular, all the existing models 
based on factorization cannot produce a large longitudinal polarization fraction, Tl/T = 0.74 it 
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NRSX model has the largest value of x (see Table VI), x — lA which deviates most from 
the value of 1.19 , and hence the smallest value oITl/V. As noted in [32], some information 
on the form factors Af^* and V^^* at = can be inferred from B K*'y decays. 

It is instructive to compare the predictions of the BSW and NRSX models for B — > 
since their B K'^*^ form factors at = are the same. Because of the dipole 
behavior of the form factors Fi, V, A2, the NRSX model predicts larger values for x, y, z and 
hence smaller values for R, Tl/T and a larger |Pp (see Table VIII). 

In short, in order to accommodate the data of S — > J/t/jK^*'> within the factorization 
framework, the form-factor models must be constructed in such a way that 

^f^K^)Mr*K^) -1-9. (5.8) 

In the hterature the predicted values of F^f-{0) spread over a large range. On the one 
hand, a large Fq^/^(0) is preferred by the abnormally large branching ratio of the charmless 
B decay B rj'K observed by CLEO [33]. On the other hand, it cannot be too large 
otherwise the SU(3)-symmetry relation FqJ^{0) — FgJ'lO) will be badly broken. There exist 
many model calculations of Fq i{0), including the lattice one, and most of them fall into the 
range of 0.20-0.33 (for a compilation of previous model calculations of Fq^^'^(O), see e.g. [34]). 
The improved upper limit on the decay mode ^ tt'^'tt", B(B^ tt'^tt") < 0.84 x 10~^ 
obtained recently by CLEO [35] implies F(fi"(0) < 0.33 or even smaller [36]. Therefore, even 
after SU(3) breaking is taken into account, it is very unlikely that F^f{0) can exceed 0.40 . 
Our best guess is that the original BSW values, F^[{0) = 0.33 and F^f{0) = 0.38 [10,11] are 
still very plausible. Taking Ff^{0) — 0.38 and using the dependence imphed by the LCSR 
(or NS, LF models), we find F^^{'m?j/^) ^ 0.70 and hence \a2\{B J/ipK) ^ 0.26 ± 0.02 
followed from Eq. (5.1). 



0.07. Various possibilities of accommodating this large Fl/F via nonfactorizable effects have been 
explored in [31,2,27]. The new CLEO [29] and CDF [30] data for Fl/F are smaller than the previous 
values. As a result, there exist some form-factor models which can explain all the three quantities 
R,Fl/F and \P\^ (see Table VIII). 
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B. Extraction of 02/01 and 02 from B D^*^Tr{p) 

The effective coefficient 02 and its sign relative to Oi can be extracted from class-Ill 
decays B~ — > D(*)'^7r~(p~) in conjunction with the class-I ones 5° — > D(*)+7r~(p~), as the 
former involve interference between external and internal VF-emission diagrams, while the 
latter proceed through the external VT-emission. Unlike the determination of Oi, there is no 
analogous differential semilcptonic distribution that can be related to the color- suppressed 
hadronic decay via factorization. Since the decay constants fn and Jd* are still unknown, 
the results for 02/01 determined from the ratios it!i^2 and it!3^4 of charged to neutral branching 
fractions [see Eq. (2.27) for the definition] are normahzed at fo — 200 MeV and fo* — 230 
MeV, respectively (Table IX). We see that 02/01 varies significantly from channel to channel 
and its value is mainly governed by Ri and R3. ^ Combining 02/01 with Table II for oi 
yields the desired results for 02 as shown in Table X. It is well known that the sign of 02 is 
positive because of the constructive interference in B~ — > D^*^^T:~{p~), which in turn implies 
that the ratios it!i, • • • , i?4 are greater than unity. 



C. Upper limit on 02 from D^tt^ 

Prom the last subsection we learn that the sign of 02/01 is fixed to be positive due to 
the constructive interference in the class-Ill modes B~ — >■ D^*^^7i~{p~), but its magnitude 
is subject to large errors. It is thus desirable to extract 02 directly from class-II modes, e.g. 

^ D(*)°7r°(p°). Although only upper limits on color-suppressed decays are available at 
present, the lowest upper limit L>°7r°) < 1.2 x 10~^ [18] can be utilized to set a 

stringent bound on 02. Neglecting 1^-exchange and final-state interactions for the moment, 
we obtain 

, 1/^ r.^r.f 0.373 \ /200MeV\ 

W(B^*)<0.29U^ ^— . (5.9) 



^The data of i^i, • • • , i?4 are taken from the Particle Data Group (PDG) [20]. Recently, CLEO 
has reported a new measurement of B ^ D*7r and obtained R3 = 1.55 ± 0.14 ± 0.15 [37], to be 
compared with i?3 = 1.67 ± 0.19 employed in Table IX. 
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The limit on 02 in various form-factor models for F^'^ is shown in Table XL 

We have argued in passing that final-state interactions (FSI) play a minor role in hadronic 
B decays, especially class-I modes. In order to have a concrete estimate of FSI, we decompose 
the physical amplitudes into their isospin amplitudes 

A{B- D\-)fsi = V3A3/2e'^3/2^ (510) 

where we have put in isospin phase shifts and assumed that inelasticity is absent or negligible 
so that the isospin phase shifts are real and the magnitude of the isospin amplitudes A1/2 and 
^43/2 is not affected by FSI. The isospin amplitudes are related to the factorizable amplitudes 
given in Eqs. (2.4), (2.8) and (2.11) by setting S1/2 = S3/2 = 0. Writing 

C = ^ W.«2(m| - ml)fnFi-{mll (5.11) 

for color-allowed and color-suppressed tree amplitudes, respectively, it is straightforward to 
show that 

A(^ ^ D\\si = A{^ ^ D\') + (e^(5i/2-53/.) _ 1) ^ 

A{^ ^ £'+7r-)Fsi = A{'^ ^ D+n-) + (^e'(5i/2-53/2) _ 1) ^ (512) 

where A(^ D^n^) = -C/v^, A(^ D+tt") = T, and we have dropped the overall 
phase e"^3/2 Taking oi = 1 and 02 = 0.25 as an illustration, we plot in Fig. 1 the effect of 
FSI on r(!B° Dtt) versus the isospin phase shift difference using the NRSX form-factor 
model. We see that FSI will suppress the decay rate of 5^ D^n~ slightly, but enhance 
that oiW ^ significantly, especially when \Si/2 — 53/2I is close to the current limit 

19° [cf. Eq. (4.2)]. This is understandable because the branching ratio of ^ D^tt^ in the 
absence of FSI is much smaller than that of — > D'^tt". Therefore, even a small amount of 
FSI via the D'^n" intermediate state will enhance the decay rate of — > D°7r° significantly. 
Fig. 2 displays the change of the upper limit of 02 in the NRSX model with respect to the 
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phase shift difference, where we have set ai — 1. Evidently, the bound on 02 becomes more 
stringent as |5i/2 — ^3/2 1 increases; we find a2{B — > Dn) < 0.29 x (200 MeV/ fo) in the absence 
of FSI and 02 < 0.21 x (200MeV//i5) at \Si/2 - ^3/2! = 19° (see Table XI for other model 
predictions) . 

D. Sign of a2{B J/^K^*)) 

Although the magnitude of 02 extracted from B J/i/jK^*^ has small errors compared 
to that determined from the interference effect in S — > D7r{p), its sign remains unknown. 
Since a2{B — > Dtt) is positive in the usual sign convention for oi, it is natural to assign 
the same sign to the J/il^K^*^ channel. It has been long advocated in [38] that the sign of 
02 (-B J/ipK) predicted by the sum rule analysis is opposite to the above expectation. 
However, we believe that a negative sign for a2{B — >■ J/ipK) is very unlikely for three main 
reasons: (i) Taking \a2{B J/i/jK) \ — 0.26 as a representative value and using cf^ — 1.149, 
cf^ — —0.325 from Eq. (2.15), we obtain two possible solutions for the nonfactorizable term 
X2iB J/i^K) [see Eq. (1.6)]: X2 = 0.18 and X2 = -0.28. Recall that X2{B Dn) is 
positive and of order 0.15 [6]. Though the energy release in B ^ J/ipK is somewhat smaller 
than that in the Dn mode, it still seems very unlikely that X2 will change the magnitude 
and in particular the sign suddenly from the Dtt channel to the J/ipK one. To make our 
point more transparent, we note that X2 has the expression: 

X2{B ^ J/i^K) = 4^^'^/'^) + ^ e?''''/'^\ (5.13) 

Cl 

where the parameters Ss and Si are defined in Eq. (1.4). Since Ci ^ 02, it is evident that 
X2 is dominated by the parameter Ss originated from color octet-octet currents; that is, the 
nonfactorized term X2 is governed by soft gluon interactions.^ Therefore, |x2| should become 
smaller when the energy released to the final-state particles becomes larger, for example. 



^In the large- A^c limit, ei is suppressed relative to eg by a factor of Nc [4]. Numerically, £i(/i) = 
-0.07 ± 0.03 and esin) = 0.13 ± 0.05 at /x = 4.6 GeV are found in [39] by extracting them from 
the data. However, it has been shown in [40] that £i(/u) is not necessarily smaller than SsilJ-), but 
this will not affect the conclusion that X2 is dominated by the eg term. 
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\X2{B 077)1 < 1X2 (-D ^ Ktt)]. It is natural to expect that \e8{B £)7r)| < \e8{B 
J/ipK)\ < \es{D KTr)\ and hence \x2iB ^ Dtt)] < \x2iB J/ipK)\ < \x2{D 
Ktt)\ as the decay particles in the latter channel are moving slower, allowing more time 
for involving soft gluon final-state interactions. Because X2(-D Kn) ~ — |, the solution 
X2{B — >■ J/ipK) = 0.28 is thus not favored by the above physical argument, (ii) Relying 
on a different approach, namely, the three-scale PQCD factorization theorem, the authors 
of [41] are able to explain the sign change of X2 from B — > Dn to D — > Kn, though the 
application of PQCD to the latter is only marginal. The same approach predicts a positive 
02 for B J/'il)K^*\ as expected [42]. (iii) The existing sum rule analysis does confirm the 
cancellation between the l/N^ Fierz term and X2 for the charmed decay D Kir [43], but 
it also shows that the cancellation persists even in hadronic two-body decays of B mesons 
[44,38,45]. For example, the light-cone QCD sum rule calculation of nonfactorizable effects 
in B^ £)°7r° in [45] yields a negative X2 and 02, which is in contradiction with experiment. 
This means that care must be taken when applying the sum rule analysis to the B decays. 
Indeed, there exist some loopholes in the conventional sum rule description of nonleptonic 
two-body decays (see also the comment made in [41]), a challenging issue we are now in 
progress for investigation. 



E. Effective Nf 

Since ci ^ C2, the effective coefficient 02 is sensitive to the nonfactorizable effects, and 
hence it is more suitable than oi for extracting Nf [strictly speaking, {Nf)2\ the effective 
number of colors defined in Eq. (2.13), or the nonfactored term X2- Although we have argued 
before that a2{B J/^/sfW) ^ 0.26 and a2{B Dn) < a2{B J/ipK^*^), it is safe to 
conclude that a2 lies in the range of 0.20-0.30. Using the renormalization scheme and scale 
independent Wilson coefficients cf^ — 1.149 and cf^ — —0.325 [cf. Eq. (2.15)], it follows that 

Nf ^ {1.8-2.2), or X2 (0.12 - 0.21), (5.14) 

recalling that X2{D Kn) ~ — |. Therefore, Nf for {V — A)(y — A) 4-quark interactions 
is of order 2 . If xi = X2, the corresponding Oi is found to be in the range of 0.97-1.01 . 
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VI. CONCLUSION 



Using the recent experimental data of S ^ D^*\n, p), B ^ D^*W^*\ and B J/tpK^*^ 
and various model calculations on form factors, we have re-analyzed the effective coefficients 
ai and 02 and their ratio. Our results are: 

• The extraction of ai and 02 from the processes B D^*^D'f^ and J/t/jK'^*^ is contam- 
inated by QCD and electroweak penguin contributions. We found that the penguin 
correction to the decay amphtude is sizable for B — > DDg, but only at the 4% level 

for B D*Ds, D^*W;, J/^K^*\ 

• The model-dependent extraction of Oi from B — > D^*^7r{p) is more reliable than that 
from B £)(*)£)(*) as the latter involve uncertainties from penguin corrections, un- 
known decay constants Jd^, Jd* and the poor precision of the measured branching 
ratios. 

• In addition to the model-dependent determination, ai has also been extracted in a 
model-independent way based on the observation that the decays B D'^*^h can 
be related by factorization to the measured semileptonic differential distribution of 
B D'^*Hd at = m^. The model- independent results ai{^ D~^7r~) = 0.93 ± 
0.10, ai(^ D+p-) = 0.95 ± 0.12 and ai(S^ ^ D*+t:-) = 1.09 ± 0.07 should be 
reliable and trustworthy. More precise measurements of the differential distribution 
are needed in order to improve the model-independent determination of a\ and to have 
a stringent test of factorization. 

• Armed with the model- independent results for ai, we have extracted heavy-to-heavy 
form factors from B D^*h: Ff^^{ml) = 0.66 ± 0.06 ± 0.05 and A^^* {ml) = 
0.56 ± 0.03 ± 0.04, where the first error is due to the measured branching ratios and 
the second one due to quark-mixing matrix elements. Form factors at other values of 
q'^ are given in Eq. (4.11). 

• Based on the assumption that 0,1 derived from B D'^*\'k, p) and from B — >■ D^*^D^/^ 
is the same, it is possible to extract the decay constants fo^ and fo* in an essentially 
model-independent way from the data. We found fo^ ~ fo* ~ C(260MeV) with large 
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errors. However, this extraction suffers from the uncertainty that we do not know how 
to estimate the violation of the above assumption. 

• By requiring that 02 extracted from J/'il)K and J/'il)K* channels be similar, as im- 
phed by the factorization hypothesis, B K^*^ form factors must respect the relation 
Ff^{mj/^) l.QAf^^rrij/^). Some existing models in which Ff^{0) is close to 
Af^*{0) and/or Ff^ does not increase with q'^ faster enough than Af^* are ruled 
out. Form factors A^'^* and V^^* can be inferred from the measurements of the 
fraction of longitudinal polarization and the P-wave component in P — > J/ipK*. 
For example, the central values of the CLEO data for these two quantities imply 

^ 1-2 and V^^*{m%^)/Af^*{m'j/^) ^ 1.5. We conjecture 
that F^^{m]/^) « 0.70 and hence \a2{B J/t/jK^*^)] ^ 0.26 ± 0.02. 

• We have determined the magnitude and the sign of 02 from class-I and class-Ill decay 
modes of S — > D^*^7r{p). Unhke 02 extracted from B — > J/'^K, its determination 
from D^*^tt{p) channels suffers from a further uncertainty due to the unknown decay 
constants fo and fo* ■ A stringent upper limit on 02 is derived from the current bound 

on B^ D^TT^ and it is sensitive to final-state interactions. We have argued that the 
sign of 02 (-B J/ipK) should be the same as a2{B Dtt) and that a2{B Dtt) ^ 
a2{B J/iPK). 

• For 02 in the range of 0.20-0.30, the effective number of colors N^^ is in the vinicity 
of 2. 
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TABLE I. Form factors at zero momentum transfer and pole masses, whenever available, in 
various form-factor models. 





BSW 


NRSX 


LF 


NS 


Yang 


LCSR 




0.690/6.7 


0.58 


0.70/7.9 


0.636 






Ff^iO)/mpoie 


0.690/6.264 


0.58 


0.70/6.59 


0.636 








0.623/6.264 


0.59 


0.73/6.73 


0.641 






Af^*{0)/mpoie 


0.651/6.73 


0.57 


0.682/7.2 


0.552 






A^^\0)/mpoie 


0.686/6.73 


0.54 


0.607/7.25 


0.441 






V^^'{0)/mpole 


0.705/6.337 


0.76 


0.783/7.43 


0.717 






F(f'^(0)/mpoie 


0.333/5.73 


0.333/5.73 


0.26/5.7 


0.257 


(see text) 


0.305 


1 \^J/ "^pole 


0.333/5.3249 


0.333/5.3248 


0.26/5.7 


0.257 


0.29/5.45 


0.305 


A^''{0)/mpoie 


0.281/5.2789 


0.281/5.2789 


0.28/5.8 


0.257 


(see text) 


0.372 


Af''{0)/mpoie 


0.283/5.37 


0.283/5.37 


0.203/5.6 


0.257 


0.12/5.45 


0.261 


A^''{0)/mpoie 


0.283/5.37 


0.283/5.37 


0.177/6.1 


0.257 


0.12/6.14 


0.223 


V'PiO) /mpole 


0.329/5.3249 


0.329/5.3248 


0.296/— 


0.257 


0.15/5.78 


0.338 


F^^{0)/mpoie 


0.379/5.3693 


0.379/5.3693 


0.34/5.83 


0.295 


(see text) 


0.341 


Ff^iO)/mpoie 


0.379/5.41 


0.379/5.41 


0.34/5.83 


0.295 


0.36/5.8 


0.341 


A^^'{0)/mpoie 


0.321/5.89 


0.321/5.89 


0.32/5.83 


0.295 


(see text) 


0.470 


Af^*iO)/mpoie 


0.328/5.90 


0.328/5.90 


0.261/5.68 


0.295 


0.18/6.1 


0.337 


A^^'{0)/mpoie 


0.331/5.90 


0.331/5.90 


0.235/6.11 


0.295 


0.17/6.04 


0.283 


V^^\0)/mpUe 


0.369/5.41 


0.369/5.41 


0.346/10.5 


0.295 


0.21/5.95 


0.458 
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TABLE 11. The effective parameter ai extracted from B —>■ (vr", p~) using different 

form-factor models. The first error comes from the experimental branching ratios shown in the last 
column and the second one from the B meson lifetimes and quark-mixing matrix elements. 







BSW 


NRSX 


LF 


NS 


Br(%) [20] 




-» D+TT- 


0.89±0.06 ± 0.07 


1.06±0.07 ± 0.08 


0.87±0.06 ± 0.07 


0.96±0.06 ± 0.07 


0.30 ± 0.04 


5° 


D+p- 


0.91±0.08 ± 0.07 


1.06±0.09 ± 0.08 


0.89±0.08 ± 0.07 


0.97±0.09 ± 0.08 


0.79 ±0.14 


S° 


D*+iT- 


0.98±0.04±0.08 


1.03±0.04±0.08 


0.83±0.03 ±0.06 


0.95±0.04 ± 0.07 


0.276 ±0.021 


S° 


D*+p- 


0.86±0.21 ± 0.07 


0.92±0.23±0.07 


0.74±0.18±0.06 


0.85±0.21±0.07 


0.67 ±0.33 


Average 


0.94±0.03 ± 0.07 


1.04±0.03 ± 0.08 


0.85±0.03 ± 0.07 


0.95±0.03 ± 0.07 





TABLE in. The effective parameter ai extracted from B D^*^Ds using different form-factor 
models. Penguin corrections to ai [see Eq. (2.18)] are included. The first error comes from the 
experimental branching ratios shown in the last column and the second one from the B meson 
lifetimes and quark- mixing matrix elements. The value of ai determined from B — > D^*^Ds and 
B D^*^D* should be multiplied by a factor of (240 MeV//£.J and (275 MeV//£.j), respectively. 







BSW 


NRSX 


LF 


NS 


Br(%) [20] 




^ D+DJ 


0.97 ±0.18 ±0.08 


1.12 ± 0.21 ± 0.09 


0.98 ±0.18 ±0.08 


1.05 ±0.20 ±0.08 


0.8 ±0.3 


B- 


- DOD- 


1.20 ±0.18 ±0.09 


1.39 ±0.21 ±0.11 


1.21 ±0.19 ±0.09 


1.29 ±0.20 ±0.10 


1.3 ±0.4 


B° 


D*+D- 


1.29 ± 0.23 ± 0.10 


1.22 ± 0.22 ± 0.09 


1.02 ±0.18 ±0.08 


1.17 ±0.21 ±0.09 


0.96 ± 0.34 


B- 


D*°D7 


1.40 ±0.29 ±0.11 


1.33 ±0.28 ±0.10 


1.01 ±0.23 ±0.09 


1.26 ±0.26 ±0.10 


1.2 ±0.5 


Avcrage 


1.17 ±0.11 ±0.08 


1.26 ± 0.11 ± 0.09 


1.08 ±0.10 ±0.08 


1.18 ±0.11 ±0.08 






D+Dl' 


1.29 ±0.22 ±0.10 


1.33 ±0.33 ±0.10 


1.15 ±0.29 ±0.09 


1.29 ±0.32 ±0.10 


1.0 ±0.5 


B- 


D°Dt- 


1.18 ±0.26 ±0.09 


1.23 ±0.27 ±0.10 


1.06 ± 0.24 ± 0.08 


1.19 ± 0.26 ± 0.09 


0.9 ±0.4 


B° 


D'+D*- 


0.91 ±0.16 ±0.17 


0.98 ±0.17 ±0.08 


0.86 ±0.15 ±0.07 


0.91 ±0.16 ±0.07 


2.0 ±0.7 


B- 


D*°Dt- 


1.09 ±0.20 ±0.08 


1.17 ±0.22 ±0.09 


1.03 ±0.19 ±0.08 


1.09 ±0.20 ±0.08 


2.7 ± 1.0 


Average 


1.05 ±0.12 ±0.08 


1.12 ±0.12 ±0.08 


0.98 ±0.11 ±0.08 


1.05 ±0.11 ±0.08 
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TABLE IV. Values of Yf' defined in Eq. (4.5) in various form-factor models. 



BSW NRSX LF NS 

1.002 1.0008 1.0009 1.001 

Y* 1.008 0.993 0.974 1.008 

Yd, 1-579 1.321 1.269 1.386 

Y?> 0.309 0.400 0.432 0.376 



TABLE V. A determination of the effective parameter ai from tfie ratio (in units of GeV^) 
defined in Eq. (4.4). The data of dB/dq^(^ D+£-i?) (in units of 10"^ GeV-^) denoted by (*) 
are explained in the text. The value of ai determined from B D^*^Ds and B D^*^D* should 
be multiplied by a factor of (240MeV//Dj and (275 MeV//D*), respectively. 





^B(^ ^ D+£-u) 


Sh 


Ol 


^B(B° ^ D*+l-u) 


SI 


ai 


ml 


0.35 ± 0.06 


0.86 ±0.19 


0.93 ±0.10 


0.237 ±0.026 


1.16 ±0.16 


1.09 ±0.07 


ml 


0.33 ± 0.06 


2.39 ±0.61 


0.95 ±0.12 


0.250 ± 0.030 


2.68 ± 1.36 


1.01 ±0.26 




0.29 ± 0.06(*) 


3.24 ± 1.07 


1.01 ±0.14 


0.483 ± 0.033 


2.09 ± 0.60 


1.22 ±0.18 


2 

Tflj-)* 

s 


0.27 ± 0.06(*) 


3.33 ± 1.34 


0.92 ±0.18 


0.507 ± 0.035 


4.32 ± 1.22 


1.05 ±0.14 



TABLE VL Form factors ^, Af^" , A^^* , V^^* and the ratio Z [see Eq. (5.5)] at q'^ = m^j 
in various form-factor models. 





BSW 


NRSX 


LF 


NS 


Yang 


LCSR 


^^^*(-^./^) 


0.56 
0.45 
0.46 
0.55 
1.08 


0.84 
0.45 
0.63 
0.82 
1.60 


0.66 
0.37 
0.43 
0.42 
1.36 


0.52 
0.39 
0.48 
0.51 
1.32 


0.50 
0.24 
0.31 
0.40 
1.04 


0.62 
0.43 
0.45 
0.86 
1.40 
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TABLE VII. The effective parameter |a2| extracted from B — > J/tpK^*^ using different 
form-factor models. Experimental branching ratios are taken from the Particle Data Group. 





BSW 


NRSX 


LF 


NS 


Yang 


LCSR 


Br(10-3) [20] 


B+ -* J/iIjK+ 


0.34±0.03 


0.23±0.02 


0.29±0.03 


0.37±0.03 


0.38±0.03 


0.30±0.03 


0.99 ±0.10 


B° -» J/ipK° 


0.33±0.03 


0.22±0.02 


0.28±0.03 


0.36±0.04 


0.37±0.04 


0.30±0.03 


0.89 ±0.12 


Average 


0.33±0.03 


0.22±0.02 


0.29±0.02 


0.36±0.03 


0.37±0.03 


0.30±0.03 




B+ J/ipK*+ 


0.20±0.02 


0.22±0.03 


0.26±0.03 


0.25±0.03 


0.40±0.05 


0.20±0.02 


1.47 ±0.27 


B° J/xpK*° 


0.20±0.02 


0.22±0.02 


0.26±0.03 


0.25±0.03 


0.40±0.04 


0.20±0.02 


1.35 ±0.18 


Average 


0.20±0.02 


0.22±0.02 


0.26±0.02 


0.25±0.02 


0.40±0.04 


0.20±0.02 





TABLE VIII. The ratio of pseudoscalar to vector meson production R, the longitudinal po- 
larization fraction Tl/T, and the P-wave component |Pp m B J/tpK^*^ decays calculated in 
various form-factor models using the factorization hypothesis. 





BSW 


NRSX 


LF 


NS 


Yang 


LCSR 


Experiment 
CLEO [29] CDF [30] 


R 


4.15 


1.58 


1.79 


3.15 


1.30 


3.40 


1.45 ± 0.26 


1.53 ±0.32 




0.57 


0.36 


0.53 


0.48 


0.42 


0.47 


0.52 ± 0.08 


0.65 ±0.11 


|P|2 


0.09 


0.24 


0.09 


0.12 


0.19 


0.23 


0.16 ±0.09 
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TABLE IX. Extraction of a^jax from B D^*)'k{p) decays in various form-factor models. 
The values of 02/01 determined from the ratios R\^2 and i?3,4 of charged to neutral branching 
fractions [see Eq. (2.27) for the definition] should be multiplied by a factor of (200MeV//D) and 

(230MeV//D*)> respectively. 





BSW 


NRSX 


LF 


NS 


Yang 


LCSR 


Expt. [20] 




O.SOibO.ll 


0.26±0.10 


0.40±0.15 


0.39±0.15 


0.36±0.16 


0.33±0.12 


1.77 ±0.29 


R2 


0.61±0.33 


0.46±0.25 


0.58±0.31 


0.52±0.32 


1.07±0.58 


0.41±0.22 


L69 ± 0.38 


Average 


0.34±0.11 


0.28±0.09 


0.43±0.13 


0.43±0.13 


0.40±0.13 


0.35±0.11 




Rs 


0.23±0.07 


0.19±0.06 


0.31±0.09 


0.28±0.08 


0.27±0.08 


0.24±0.07 


1.67 ±0.19 


R4 


0.55±0.45 


0.64±0.52 


0.85±0.70 


0.74±0.61 


1.47±1.20 


0.61±0.50 


2.31 ±1.23 


Average 


0.24dr0.07 


0.19±0.06 


0.32±0.09 


0.29±0.08 


0.28±0.08 


0.25±0.07 





TABLE X. The effective coefficient 02 extracted from the analyses of i?j and 

— > D(*)+7r~(p~). The values determined from Ri^2 and i?3,4 should be multiplied by a fac- 
tor of (200MeV//£i) and (230MeV//D*), respectively. 

BSW NRSX LF NS Yang LCSR 

i?i & 5° ^ D+TT- 0.27±0.10 0.27±0.10 0.35±0.13 0.38±0.14 0.35±0.13 0.32±0.12 

i?2 & S° ^ D+p- 0.55±0.31 0.49±0.27 0.51±0.28 0.58±0.32 1.04±0.57 0.40±0.22 

Average 0.30±0.10 0.30±0.10 0.38±0.12 0.41±0.13 0.38±0.13 0.33±0.11 

R^kB^^ D*+Tr- 0.22±0.07 0.19±0.06 0.26±0.08 0.27±0.08 0.26±0.08 0.23±0.07 

i?4 & 5° ^ D*+p- 0.47±0.41 0.59±0.50 0.63±0.54 0.63±0.54 1.24±1.07 0.52±0.44 

Average 0.23±0.07 0.20±0.06 0.26±0.08 0.28±0.08 0.26±0.08 0.24±0.07 
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TABLE XL The upper limit on the effective coefficient 02 [multiphed by (200 MeY/fo)] inferred 
from the decay D^tt^ in the absence and presence of final-state interactions characterized by 

the isospin phase shift difference A = \di/2 ~ <^3/2l-B-+D7r- 







BSW 


NRSX 


LF 


NS 


Yang 


LCSR 


02 (with A = 


0°) 


0.29 


0.29 


0.38 


0.41 


0.38 


0.34 


02 (with A = 


19°) 


0.17 


0.21 


0.21 


0.26 


0.24 


0.22 




FIG. 1. The ratio of r{B Dtt) in the presence of final-state interactions (FSI) to that 
without FSI versus the isospin phase-shift difference. The calculation is done in the NRSX model 
[12]. 
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1^1/2-^3/21 



FIG. 2. The upper bound of the effective coefficient 02 multiphed by (200 MeV//^)) derived 
from the current hmit on D^ir^ using the NRSX model [12] versus the isospin phase-shift 

difference. 
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